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Abstract 

The density of states and other thermodynamic functions of hehum nanodroplets are calcu- 
lated for a microcanonical ensemble with both energy and total angular momentum treated as 
conserved quantum numbers. These functions allow angular momentum conserving evaporative 
cooling simulations. As part of this project, a recursion relationship is derived for the reduction 
to irreducible representations of the n'th symmetric power of the irreducible representations of the 
rotation group. These give the distribution of total angular momentum states generated by putting 
multiple quanta into a ripplon or phonon mode of the droplet, each of which is characterized by a 
angular momentum quantum number. 
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INTRODUCTION 



Spectroscopy of atoms and molecules dissolved in superfluid helium nanodroplets has 
become a field of intense activity. Such nanodroplets rapidly cool by evaporation, reaching 
final temperatures of ~ 0.37 KJi|. There is a long history of their formation and study 
mass spectroscopic methods j2,|3|. Vibrational , electronic 0,0,0], and rotational 



10 



ll| transitions have been observed. Part of the interest in this field is that helium droplets 



provide a nearly ideal matrix [12l] for isolation of highly unstable compounds. Many novel 
compounds and clusters have been synthesized and spectroscopically studied in this unique 



environment 
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15|,ll6|,ll7|. 



There have been several previously reported calculations of the evaporative cooling of 
helium nanodroplets 19], and these have predicted a terminal temperature in excellent 
agreement with that later observed experimentally jl|. These calculations, however, have 
not considered angular momentum conservation, which obviously imposes constraints on 
droplet cooling in a high vacuum environment. In order to address the question of the 
possible trapping of angular momentum in a cooling droplet, we have undertaken a new 
study of the evaporative cooling of the nanodroplets, but including angular momentum con- 
servation, using methods analogous to "Phase Space Theory" calculations of unimolecular 
dissociation |2(|. Necessary inputs to such calculations are the density of states and inte- 
grated density of states of helium nanodroplets, over the energy and angular momentum 
range sampled in the evaporative cooling trajectories of nanodroplets. 

The angular momentum resolved density of states of helium nanodroplets is not available 
in the literature. While the calculation of this quantity is largely based upon straightforward 
extensions of standard convolution methods, there is one particular point that was not. The 
distribution of angular momentum states produced by excitation of an arbitrary number of 
quanta in a (2L + 1) degenerate vibrational state of angular momentum quantum number L 
was required. This distribution is provided by the reduction of the ra'th symmetric product 
of the L'th irreducible representation of the group of transformations of a sphere, K. While 



the equivalent reduction for common molecular point groups is well known |2l| , the present 
author could not locate the general result for K, and has derived a recursion relationship 
that provides the needed reduction. This is presented in the present paper, along with a 
derived asymptotic expression, which is of the same form as the thermal distribution of a 



2 



rigid spherical top. 

The density of both ripplon and phonon modes are considered separately. In each case, 
the simple scaling of the excitation spectrum with droplet size (and thus number of helium 
atoms), allows the density of states to be derived as a unique function of a reduced energy. 
These functions have been fit to simple analytical expressions which provide excellent ap- 
proximations over the range of energy and angular momentum of interest. In particular, it 
is found that the distribution of angular momentum quantum states is of the form of the 
thermal distribution of rotation for a rigid spherical top, with a reduced energy dependent 
effective inverse "temperature" that changes slowly with energy. 



RIPPLON EXCITATIONS 



The lowest energy excitations of a pure droplet are ripplons, which are quantized capillary 
waves on the surface of the droplets. The properties of such waves on spherical drops are 
thoroughly described in the classic text on hydrodynamics by Lamb [2^. Classically, each 
ripplon mode involves a modulation of the surface of the droplet, with 

r{e, (j)) = ro+J2Re {AlmYlm{0, <^)e''^^*) (1) 

L,M 

where L > 2 and 



u;l = jL{L-l){L + 2)uJo (2) 



and a is the helium surface tension (approximated by the bulk zero temperature value of 
0.363 mJ/m^), m is the atomic mass of helium, p is the number density of helium (approx- 
imated by the bulk zero temperature value of 0.0218 A'^), and R = = 2.22N^/^ k 
with N equal to the number of helium atoms in the droplet. Given that the spectrum of 
ripplon modes is proportional to a common factor that contains the size dependence, the 
thermodynamic functions of the ripplon modes are universal functions of a reduced energy 
which has unit value equal to = Kluq = 3.77/-\/iV k^K. 

The integrated total state count, Ne, and the total density of states, pe, can be calcu- 
lated from the spectrum of ripplon mode frequencies using the Beyer-Swinehart direct count 
method [2^. In addition to the spectrum, one needs the degeneracy of a state with ul 
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quanta in the ripplon mode with quantum number L. This is equal to the number of distinct 
ways of putting ul identical objects into 2L + 1 bins and is equal to (/(n/,, L) = ''^'^^ji^^l' 
Define Ne{E, L) to be the integrated density of states calculated using only ripplon modes 
with angular momentum quantum numbers less than or equal to L. We can recursively 
calculate the total integrated density of states by using: 

Ne{E, L = 2) = Y^^E- hu2n2)gin2, 2) (4) 
Ne{E, L) = Y. N{E - huLUL, L - l)g{nL, L) (5) 

where the unit step function is defined by <l>(x) = 1(0) for x greater or equal to (or less than) 
zero. In the Beyer- Swinehart algorithm, the above are calculated using an energy bin size, 
5E sufficiently small that each ujl/5E can be approximated by an integer. A step size of 0.1 
reduced energy unit was used in this work. The calculated values of A''^; and pE are plotted 
on figure [T] We have fitted the calculated integrated density of states to the expression 

Ne = exp (aE^/^ + (6) 

Values of a = 2.5118 and b= -3.4098 minimized the root mean squared fractional errors 

over the interval E = 50 2500. Changing the step size in the calculation and refitting 

the integrated density of states gave values of a = 2.5118 and b = -3.4110. The leading 
power of 4/7 can be derived by equating the high temperature limit of the microcanonical 
and canonical ensemble values of energy and entropy (which is proportional to ln{pE) in the 
microcanonical case). The power of 1/7 for the correction term was selected based upon the 
slope, on a log-log plot, of NeE'^/"^ vs. E. The expression for Ne given in Eq.lHl agrees with 
the Beyer- Swinehart direct count to within ±5% while Ne increases by more than 83 order 



of magnitude over the reduced energy range 50 — 2500. Using standard expressions 24| for 



the thermodynamic functions of a microcanonical ensemble, we can derive the following: 

(7) 



Pe{E) = f^j = Ne{E) 



^aE'^/' + -bE-^/' 
7 7 



Se{E) = kj,ln{pE{E)) = 



T (F) = = ^ 28aE^°/^ + 7E e^i J_1j,3/7 .q^ 

'^^ ' [dE k^lQa^E + MhE^n -l2aE^n + h^EyT -Qh k^Aa ^ ' 



V 



'dry' 

dE 
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6 (32a3Eio/7 + 32a%E - dGa^E^/^ + lOab^E'^/'^ - iSabE^^ + b^E^^ - 7b^) 

Se, Te, and Cy are the microcannonical values of the entropy, temperature, and heat capac- 
ity due to the droplet modes as a function of reduced rippon energy. The density of states, 
Pe, can be compared with O.Slli?"^/'' exp(2.49i?^/^) derived by Brink and Stringari [l8| 
through inverse Laplace Transformation of the high temperature canonical partition func- 
tion using the stationary phase approximation. (Their reported results have been converted 
to the normalized roton energy units used in this work) . It can be seen that the exponential 
dependence of the density of states is almost identical, but there is a slight difference in the 
power of the energy dependence of the exponential prefactor. 

For isolated helium nanodroplets, angular momentum as well as energy is a conserved 
quantity, and thus the proper statistical ensemble is one that sums only over states of the 
same energy and total angular momentum. As a first step in the calculation of the properties 
of this ensemble, we have derived a regression expression that determines Nj{J, ni, L) which 
is the number of states of total angular momentum quantum number J that are generated 
from placing ni quanta in a ripplon mode with angular momentum L. We "count" such that 
each "state" is an irreducible representation with 2 J + 1 values of the projection quantum 
number. This recursion expression and its derivation is given in the Appendix to this paper. 

The individual values of Nj{J, n^, L) are quite irregular and do not appear to derive from 
any simple expression. However, the sum of states, of course, gives the simple expression: 

5:(2 J + l)Nj{J, n,, L) = g{n,, L) = ^^^r^ (11) 

Also, the mean value of the squared angular momentum has been empirically found to be 
given by a simple expression: 



< J(J+ 1) >ni,L= riL 



L{L + l) + h{nL-l] 



(12) 



If we ignored the Bose symmetry (i.e. treated the angular momenta of different quanta in 
the same Ripplon mode as uncorrelated) , we would not have the last term on the right hand 
side. The Bose symmetry increases the mean squared angular momentum because ripplon 
excitations have a greater probability of being in the same direction than if they where 
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distinguishable. This closed form expression can be used to sum a thermal distribution of 
ripplons as a function of inverse temperature in reduced energy units, jS = Ey/ k-oT. 



<J{J + 1)>=Y,<J(J+ 1) >„^,i g{nL, L) exp ( -PriL^LiL - 1)(L + 2] 

nL,L 



^L(L + 1) 

L 



2L + 1 



1 



21 L/3y/LiL~l)iL+2) _ 2 



(13) 



If we neglect the correlations induced by the Bose symmetry of the ripplon modes, we get 
the same result as above but with the last term on the right hand side omitted. 

In a previous paper, the present author used the uncorrelated expression to calculate the 
root mean squared averaged angular ripplon angular momentum as a function of droplet 
size j|2^. Figure El shows a plot of the thermally averaged value of calculated both ne- 
glecting the Bose symmetry and by explicit thermal average over the present Microcanonical 
results. It is seen that the exact calculation is modestly higher than for the uncorrelated 
prediction, but by a nearly constant factor of approximately 20%. 

Given Nj{J,ni, L), the J restricted state count and density of states can be computed 
using a modification of the Beyer-Swinehart direct count. Define Nej{E, J, L) to the in- 
tegrated density of states with total angular momentum J, but including excitation only 
in ripplon modes up to angular momentum quantum number L. The following recursion 
relation is easily derived by using the triangle rule: 

Nej{E, J, L = 2) = ^ $(£; - hu2n,)Nj{J, n^, 2) (14) 

n2 

J+J' 

Nej{E,J,L) = Yl E NEj{E-huLnL,J",L-l)Nj{J',nj,L) (15) 

nL,J' J" = \J-J'\ 

for each E, we need to iterate on L until Hujl > E to get the complete integrated density of 
states, Nej{E, J), with total angular momentum J. We have explicitly calculated Nej{E) 
for E < 200. by this procedure. It was found that for each E > 50, Nej{E, J) has a 
distribution that accurately fits the functional form: 



Nel{E, J) = iV^(E)(2J+ 1) J^^exp (-/5i(E)(J + 1/2)2) (^g) 



TX 

which is the form of the thermal distribution for a spherical rigid rotor. Figure |21 shows plots 
of the calculated distributions compared to the fitted "thermal" form for several values of 
the reduced energy. 



Assuming this "thermal" distribution of states with different rotational quantum num- 
bers, The value of (3l{E) can be calculated if we know the mean value of total angular mo- 
mentum, < J(J + 1) >, over all states up to a given energy by using < J(J + 1) >= \l3j} . 
Using the relationship given in Eq. El we have derived a recursion relationship to calculate 
this quantity without having to explicitly enumerate all the states as a function of J . Let 
A{E) be the average of the total angular momentum quantum numbers for all states up to 
reduced energy E, and A{E, L) the same quantity, but calculated only using ripplon modes 
with quantum numbers up to L. Since angular momentum in different modes are uncorre- 
lated, the total angular momentum is just the sum of angular momentum in each individual 
mode. It is straightforward to derive that 



Ne{E,L)A{E,L) = Y.g{nL,L)NEiE-hiULnL,L-l) 



X 



■n-L 



{A{E-huLnL,L-l] 



L(L + l) + ^LK-i: 



(17) 



The inverse "rotational temperature" , Pl{E), can be calculated by using the standard spher- 
ical top relation A{E) = 3/{2Pl{E)). Using this procedure, Pl{E) has been calculated for 
E = 10 — 2500.. The values of Pl{E) are plotted in figure IH along with a fit through the 
points to a power law in energy: 

Pl{E)=cE~^'' + dE-^^" (18) 

with c = 0.8680 and d = 0.9639. A fit the the values calculated using an energy bin of 0.01 
reduced units gave c = 0.8679 and 0.9759. The Pl values from this fit are also displayed on 
figure HJ but cannot be distinguished on the scale of the plot. The Pl{E) calculated with 
Eq. ITHl agrees with the calculated values to within 0.1% for E = 100 — 2500 and to within 
1% for E = 50 — 100. The reduction in accuracy in this domain is because the calculated 
values of Pl{E) oscillate around the value predicted by Eq. HSl The leading power of -8/7 
in (3l{E) can be derived by equating the high temperature limit for the mean value of the 
total angular momentum squared for the canonical and microcanonical ensemble. 

Using the standard relationships, we can use calculate the thermodynamic quantities as 
a function of reduced ripplon energy and total angular momentum of the droplet as: 



PeAE,J) = (^] =Nej 



dE 



J 



4 1 «/7 / 3 / 1\^\ d(3 



7 7 \2pL V 2J I \dE 



(19) 
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Sej{E,J) = kbln{pEj) = h 



,E^/7 + 5^1/7 + 1^(2 J + 1) + i In (f^) _ J + 



+ ln 



dE 



(20) 



1 



kbTEj{E, J) h 



dE 



7 7 



2/?i 



J 



3 /cf^V 
Wl\dE ) 



3 
2/3i 



(^) 



(21) 



If we neglect ^Ir (^)^ in Eq. EH (which ^ 3.2^-^/^ times ^aE'^^/'^ that it is added 



to), then the equation for Tgj can be shown to reduce in the high energy hmit to: 



Tej{E, J) = Te{E) + 0A8E-'/' f-|- 



J 



o 



J + 



2\ 2 



(22) 



2Pl V 2J J \2(3l 

Tej decreases monotonically with J for fixed E. Tej{E,0) = Te{E) + while for 

the RMS value of J, Tej{E, J) = Te{E). However, for J > 1.2^6/7, Tej -> IME^^/'^J-^ 
values much below Te{E) for the same value of E. 



PHONON EXCITATIONS 



Starting at higher energy than the surface ripplon modes are phonon (compressional) 
excitations of the helium droplets. These normal modes are characterized by two quantum 
numbers, n, L, where n = 1,2,... is the number of radial nodes and L is the angular 
momentum quantum number. Each normal node can be characterized by a wavenumber 
kn^L which is determined by kn^i = fn^ilR where t^^l is the n'th root of the spherical bessel 
function ji [26]. For k^^i 1 (i.e. much smaller than the wavenumber of a roton), the 
excitation angular frequency of each mode is given by uJn,L = ukn^L, where u = 236 m/ s is the 
speed of sound in helium 26|. We define a reduced energy Ep with an energy unit equal to 
the excitation energy of the lowest (n = 1, L = 0) phonon which is hun/R = /i-534A^^^/'^ GHz 
= 25.5N-^/^KkB. 

We have calculated the density of states using similar methods as for the ripplon density 
of states. The integrated density of states was fit to the functional form 



,^3/4 ^ ^^1/4 ^ f 



(23) 
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The values a = 3.3306, b = —3.5941, and / = 1.7786 reproduces the recursively calculated 
values of Ne to an accuracy of better than ±0.09% while Ne varies by more than 130 orders 
of magnitude over the reduced energy interval Ep = 25 — 500. 

We have assumed that the distribution of total angular momentum states will again be 
well approximated by the spherical top thermal distribution expression ITHl Values of I3l{E) 
have been calculated using the recursion relationship for A{E)^ Eq. El These were found 
to be well approximated by the expression 

Pl{E) = cE-''/^ + dE-'/^ (24) 

Values of c = 0.2545 and d = 0.2929 reproduce the calculated values of I3l{E) to a fractional 
accuracy of better than 0.2% over the reduced energy interval E-^ = 25 — 500. Given these 
expressions for Ne{E) and (3l{E), expressions for the various microcanonical thermodynamic 
quantities can be derived, both with and without angular momentum constraints, analogous 
the expressions given above for the ensemble of ripplon modes. 

The same procedure can be used to calculate the density of states without the restriction 
kn,L ^ 1 A, if the values of z/„ ^ are calculated using the elementary excitation energy of 
liquid helium for wavenumber /c„,l. However, this no longer results in a density of states 
that is a function only of a reduced energy, i?p, and thus would have to explicitly calculated 
for each droplet size. This has not been pursued in the present work. 

CONCLUSIONS 

In this work, we have presented simple numerical procedures that allows the calculation 
of the density of states of both ripplon and phonon excitations of a helium nanodroplet as 
a function of both reduced energy and total angular momentum quantum numbers. It is 
found that these can well approximated by simple analytical expressions. It has been found 
that the distribution of rotational total angular momentum, for each energy, closely follows 
that of rigid spherical top. The microcanonical expressions for the ripplon density of states 
has been used by the author and collaborator for statistical evaporative cooling calculations 
of both pure and doped helium nanodroplets, conserving angular momentum. That work 
will be presented in a future publication. 
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APPENDIX: REDUCTION OF SYMMETRIC PRODUCTS OF THE ROTATION 
GROUP 

In order to compute the density of ripplon states as a function of both energy and total 
angular momentum quantum number, one needs to determine the distribution of states 
of different angular momentum generated by quanta in each ripplon mode with mode 
angular momentum quantum number L. We use the fact that a set of states with total 
angular moment quantum number L and with projection M = —L, —L + 1 ... L is an 
irreducible representation of the rotation group of the sphere, K (or the double group of 
rotations is half integer angular momentum is allowed for). Below, when we speak of a 
"state" with a given total angular momentum quantum number, J, it will be implicit we 
mean a degenerate set of (2 J + 1) eigenstates with the allowed range of projection quantum 
numbers. For the states generated by multiple excitations in different modes, the total 
angular momentum distribution is found by reduction of the direct product representation 
of the distributions in each mode. The reduction of a direct product representation for the 
sphere gives the well known triangle rule, i.e. products of states with quantum numbers Ji 
and J2 give one state each with J = | Ji — J2I, l-'^i — -^2! + 1, • • • , -^1 + -^2- The triangle rule 
can be applied recursively to determine the total number of states with each total J once 
we know the distribution of number of states with each J quantum number for excitation 
of riL quanta in a harmonic ripplon mode with angular momentum quantum number L. 

The distribution of states produced by multiple excitation in a degenerate mode is not 
give by the direct product of that mode with itself rii times, but by the states given from 
the reduction of the symmetric ul power product of that mode. This is the meaning of the 
widely cited statement that vibrational quanta are "Bosons" . For the point groups relevant 
to the spectra of rigid molecules, the reduction of such symmetric products can be found in 
a standard text on vibrational spectroscopy Q- However, the author has not been able to 
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locate equivalent expressions for the reduction of general symmetric powers of the rotation 
group. In this appendix, we give a recursion relation that has been used to calculate these 
up to high values of L and til. 

Each state produced by the symmetric direct product of quanta in mode with total 
angular momentum quantum number L can be represented by a set of mj quantum numbers 
(integer or half integer depending upon L) such that 

— L < mi < 1712 • • • < ^riL ^ L (A. 25) 

We will calculate a recursive relationship for NM{M,nL, L) which is the number of n^'th 
power symmetric product states of mode L such that the sum of the quantum numbers mj 
equals M. For = 1, Nm{M, 1,L) = 1 if \M\ < L and zero otherwise. Note that if we 
restrict = L, then the total number of such states that contribute to NM{M,nL, L) is 
equal to Nm{M—L, ul — I, L). If = L— 1, then because the rest of the quantum numbers 
are restricted to only 2L values, the number of states that contribute to NmIMjUl, L) is 
equal to Nm{M — (L — 1) + {n — l)/2, rii — 1, L — 1/2). The last term in the first argument 
arises from the "shift" in the projection quantum number when mapping the sum over the 
2nL values to the symmetric direct product for a mode with angular momentum quantum 
number L — 1/2. For the general case of = L — k, the number of such states contributing 
to Nm{M, m, L) is equal to Nm{M - {L - k) + k{nL - 1) /2, - 1, L - k/2). Taking these 
contributions into account as well as that < k < 2L and that NM{M,nL, L) = if 
\M\ > ulL we can write: 

^max 

NM{M,nL,L)= NMiM-L + k{nL + l)/2,nL-l,L-k/2) (A.26) 

Where kmin is the maximum value of and —M + 2L — ulL and k^ax is the minimum value 
of 2L and L — [M/til). We can then calculate Nj{J, n^, L), the number of states with total 
angular momentum quantum number J by Nj {J,nL,L) = Nm {Jj^l, L) — Nm {J +^,nL,L), 
as is done in many introductory Quantum texts when calculating the terms produced by 
a given atomic configuration. It is apparent, since all terms in Eq. IA.26I shift in their 
first argument by one when the first argument on the left shifted by one that exactly the 
same recursion relationship can be applied to directly calculate the Nj{J,nL, L) values. 
However, in such a calculation, one must replaces Nj{—J,nL,L) by —Nj{J~ l,ni,L) and 
Nj{— 1/2, riL, L) = when these negative argument values arise in the recursion. Also, 
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the recursion is started by using Nj{J, 1,L) = S{J,L). It is also easily demonstrated that 
Nj{J,2,L) = (1/2) (l + (-l)-^+2^) for J < 2L (the factor of 2L is included to make the 
result correct even in the case of half integer quantum numbers). Note that it is possible 
to calculate each successive value of 11^ by replacement in a single stored two dimensional 
matrix by looping down from the highest L value. 

The calculation of the Nj{J,nL, L) values becomes quite memory demanding for large 
values of nL,L. This is in part because they must be generated with increasing values of rii, 
but must be used in the calculation of the density of states in order of increasing L, thus 
requiring a three dimensional matrix with ^ |(^max-Z^max)^ elements, where nmax and 
Lm&x are the highest values of and L needed. However, it has been found that for large 
hl and L values the distribution of Nj{J, nL,L) approaches that of a rotational distribution 
of a spherical top, i.e. 



Nj{J, UL, L) ^ g{nL, L){2J + 1)^^^!^^ exp {-(3{nL, L){J + 1/2)^) (A.27) 

where g{nL, L) = ^('^^ji^'^l' is the total number of states (including the 2J+1 degeneracies) 
generated with J = . . . n^L. This functional form can be justified when it is recognized 
that the central limit theorem j23| suggests that A^a/(M, ni,L) should approach a Gaussian 
distribution in M for large rii values and that Nj{J, riL.L) is the derivative of Nm{J, til.L) 
in that limit. We can determine the dependence of (3{nL,L) by equating the "thermal 
average" value of J{J + 1) to the exact expression, Eq. which gives 

PiriL, L) = . (A.28) 

2nL[L{L + l) + \L{nL-l)\ 

For rti = 14 and L > 10, the maximum error in this approximation to Nj{J,nL, L) is less 
than 5% of the peak value for each L. 
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FIG. 1: Plot of the natural log of the integrated density of states and density of states in reduced 
energy units. The integrated density is slightly higher than the density of states. 
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FIG. 2: Plot of the Root Mean Squared total thermal angular momentum in ripplons as a function 
of reduced temperature. It can bee seen that the effect of the Bose correlations of the ripplon 
excitations increases the mean angular momentum by about 20% over the uncorrelated result. 
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FIG. 3: Plot of normalized integrated density of states as function of total angular momentum 
quantum number for three different values of the reduced energy. Also plotted are the distributions 
predicted by fits of these distributions to a thermal distribution of a spherical top, as given by 

Eq.nni 
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FIG. 4: Plot of I3l{E) as function of reduced energy. The values of /?l where calculated both 
by from the mean value of the squared total angular momentum, as calculated by the recursion 
relations, and using Eq. 1181 These two values agree to better than the resolution of the plot. 
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